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We study Polyakov loop correlations and spatial Wilson loop at finite temperature in two-flavor 
QCD simulations with the RG-improved gluon action and the clover-improved Wilson quark action 
on a 16"^ X 4 lattice. From the line of constant physics at mps/my ~ 0.65 and 0.80, we extract 
the heavy-quark free energies, the effective running coupling gcff (r) and the Debye screening mass 
\ moiT) for various color channels of heavy quark-quark and quark-anti-quark pairs above the critical 

. temperature. The free energies are well approximated by the screened Coulomb form with the 

appropriate Casimir factors at high temperature. The magnitude and the temperature dependence 
I of the Debye mass are compared to those of the next-to-leading order thermal perturbation theory 

. and to a phenomenological formula in terms of g^siT). We make a comparison between our results 

with the Wilson quark action and the previous results with the staggered quark action. The spatial 
fT^ ' string tension is also studied in the high temperature phase and is compared to the next-to-next- 

leading order prediction in an effective theory with dimensional reduction. 
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O ■ I. INTRODUCTION 

o 

Recent relativistic heavy-ion experiments have revealed various remarkable properties of QCD at finite temperatures 
• and densities, suggesting the realization of the QCD phase transition from the hadronic matter to the quark-gluon 
' plasma (QGP) In order to extract unambiguous signals for the transition from the heavy-ion experiments, it is 
indispensable to make quantitative calculation of the thermal properties of QGP from first principles. Currently, the 
lattice QCD simulation is the only systematic method to do so. By now, most of the lattice QCD studies at finite 
' , temperature and chemical potential have been performed using staggered quark actions with the fourth-root trick 
of the quark determinant, which require less computational costs than others. However the lattice artifacts of the 
(-H [ staggered quark actions are not fully understood. Therefore, it is important to compare the results from other lattice 
!• • quarks such as Wilson quark actions to control and estimate the lattice discretization errors. 
. 5^ I Such a study at finite temperature (T ^ 0) and zero chemical potential (/x, = 0) has been initiated several years 
^ . ago using the Iwasaki (RG) improved gauge action and the Nf = 2 clover improved Wilson quark action by the 
■ CP-PACS Collaboration The phase structure, the transition temperature and the equation of state have been 

investigated in detail, and also the crossover scaling around the chiral phase transition has been tested. In contrast 
to the case of staggered quark actions, the subtracted chiral condensate in the standard Wilson quark action [3] and 
in the clover-improved Wilson quark action [5] shows the scaling behavior with the critical exponents and scaling 
function of the three-dimensional 0(4) spin model. This suggests that the lattice QCD with Wilson- type quarks is 
in the same universality class as the 0(4) spin model, as expected from the effective sigma model analysis [5!, 
Moreover, extensive calculations of various physical quantities at T = such as the light hadron masses have been 
carried out using the same action 0, § . 

Since a lot of experimental results are obtained by the heavy-ion collisions and numbers of technical progresses in 
treating system at finite baryon density on the lattice have been made after the studies by the CP-PACS Collaboration, 
it is worth while to revisit the QCD thermodynamics with Wilson quark actions. In particular, it is essential to perform 
simulations along the lines of constant physics (LCP) to clearly extract the temperature- and density-dependences. 
As a first step in this direction, we carry out simulations of Nf = 2 QCD on an x Nt — 16'^ x 4 lattice at 
mps/my — 0.65 and 0.80 in the range T/Tpc ^ 0.76-4.0, where mps (my) is the pseudo-scalar (vector) meson mass 
and Tpc is the pseudo-critical point along the LCP. Among various topics which can be studied using the above 
configurations, we will focus on two subjects in this paper: the free energy between heavy quarks and the spatial 
string tension. They are the fundamental quantities to characterize the perturbative and non-perturbative properties 
of the hot QCD medium. 

The heavy-quark free energy was recently studied by lattice simulations in the quenched approximation [l^, [ll[ , 
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and in full QCD with the staggered quark action [T^, [T^ and with the Wilson quark action • An analytic study 
based on the thermal perturbation theory was also reported p^ . It was pointed out that the interaction between 
heavy quarks is intimately related to the fate of charmoniums and bottomoniums in QGP created in relativistic 
heavy- ion collisions [l6| . In this paper, we make a systematic study of the free energy between a quark (Q) and 
an antiquark (Q) in the color singlet and octet channels, and between Q and Q in the color anti-triplet and sextet 
channels. We adopt the Coulomb gauge fixing for the gauge-non-singlet free energies. By fitting the numerical results 
with a screened Coulomb form, we extract an effective running coupling and the Debye screening mass in each channel 
as functions of temperature. We also study the "force" between heavy quarks defined through the derivative of the 
free energy with respect to the inter-quark distance. We show that (i) the free energies in the different channels at 
high temperature (T > 2Tpc) can be well described by channel-dependent Casimir factors together with the channel- 
independent running coupling gcf{{T) and Debye mass moiT), (ii) the next-to- leading order result of the Debye mass 
in thermal perturbation theory agrees better with the lattice m]j{T) data than that of the leading order, (iii) gcf[{T) 
and moiT) satisfy the leading order formula, moiT) = ^1 + Nf /6 geff{T) T for T > l.bTpc, so that most higher order 
and non-perturbative effects on the Debye mass is likely to be absorbed in gcs{T), and (iv) there is a quantitative 
discrepancy in tod (T) between our results using the Wilson quark action and those using the staggered quark action 
[l^ [isj even at T ~ 4Tpc, and (v) the Casimir scaling law valid above Tpc is violated below Tpc, in particular in the 
color octet channel. 

We also extract the spatial string tension a-^ at T > Tpc from the spatial Wilson loop. We show that ^ctT has 
approximate linear increase with T, as previously reported in quenched studies [13, fisl. Il9l. [20j . We find that our 
result for cTs agree quantitatively well with the analytic result based on a dimensionally reduced effective theory [2l[ . 

This paper is organized as follows: In Section |TT1 we present our lattice action and simulation parameters, and 
discuss the line of constant physics. Results of numerical simulations for the heavy quark free energies are shown 
in Section IIIII The effective running coupling and Debye screening mass are extracted from the heavy quark free 
energies and are compared with the analytic results of thermal perturbation theory and with the numerical results 
with staggered quark action. The force between heavy quarks are also discussed in this section. In section IIV| we 
show the spatial string tension obtained from the spatial Wilson loop and its comparison to analytic results. The 
paper is concluded in Section |Vl We discuss the fit range dependence of the free energies in Appendix]^ and tabulate 
numerical data of the free energies in Appendix IB] 



II. SIMULATIONS WITH A WILSON-TYPE QUARK ACTION 



A. Lattice action 

We employ the RG-improved gauge action ^251 and the Nf = 2 clover-improved Wilson quark action [2^ defined 

by 

S = Sg+Sg, (1) 

/=1,2 x,y 

where /3 = 6/5^, c\ = —0.331, cq — 1 ~ 8ci and 

Dx,y = 6xy - K'Y{{1 - Jfj_)Ux,f,Sx + f,,y + {I + Jfj.)U^^^Sx,y+fj.} - SxyCSwK^^Crfj^^Ffj,^. (4) 

Here K is the hopping parameter and F^j, is the lattice field strength, F^^ = ~ //Iiy); with /^j^ the standard 

clover-shaped combination of gauge links. For the clover coefficient csw, we adopt a mean field value using W^^^ 
which was calculated in the one-loop perturbation theory j22| . 

csw = (VF^^i)-^/^ - (1 - 0.8412/3-1)-^/^ 

The phase diagram of this action in the (/3, K) plane has been obtained by the CP-PACS Collaboration 0, 0] as 
shown in Fig. [TJ The solid line Kc{T = 0) with filled circles is the chiral limit where pseudo-scalar mass vanishes at 
zero temperature. Above the Kc{T = 0) line, the parity-flavor symmetry is spontaneously broken [23, [2^, [2^. At 
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FIG. 1: Phase diagram for RG improved gauge action and clover improved Wilson quark action for Nt = 4. 
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FIG. 2: Pseudo-scalar meson mass squared (left) and vector meson mass (right) as a function ot 1/K for several values of /3 at 
T = 0. 



finite temperatures, the cusp of the parity-broken phase retracts from the large (3 limit to a finite (3 [27ll28l|. The solid 
line Kc{T > 0) connecting open symbols represents the boundary of the parity-broken phase. The region below Kc 
corresponds to the two-flavor QCD with finite quark mass. We perform simulations in this region after investigating 
the relation between the simulation parameters (/3, K) and the physical parameters, e.g. quark mass and lattice 
spacing. 

The dashed line Kt with filled diamond represents the finite temperature pseudo-critical line determined from the 
peak of Polyakov loop susceptibility. This line separates the hot phase (the quark-gluon plasma phase) and the cold 
phase (the hadron phase). The crossing of the Kt and the Kc{T — 0) lines is the chiral phase transition point. 

B. Determination of lines of constant physics and simulation parameters 

For phenomenological applications, we need to investigate the temperature dependence of thermodynamic observ- 
ables on a line of constant physics (LCP), which we determine by mpg/my (the ratio of pseudo-scalar and vector 
meson masses at T — 0). For our purpose, we need LCP in a wider range of parameters than that studied in 
Ref. [sj. Therefore, we re-analyze the data for mpsa and my a at zero temperature in a wider range shown in Fig. [21 
where a is the lattice spacing 

0, i, \l, 3]- The thin solid lines in Fig. [3] shows our results for LCP corresponding to 
™ps/™v = 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 and 0.95. The bold solid line denoted as Kc represents the critical line, 
i.e. mps/my = 0. Our LCP's are consistent with those of l3j] in the range of the previous study. 

We also re-analyze the lines of constant T/Tpc. The temperature T is estimated by the zero-temperature vector 
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FIG. 3: Lines of constant mps/mv (solid lines) in the {l3,K) plane for mps/mv = 0.65, 0.70, 0.75, 0.80, 0.85, 0.90 and 0.95. 
Dashed lines represent lines of constant T/Tpc on A^t = 4 lattices, where Tpc is the pseudo-critical temperature corresponding 
to Kt{Nt — 4) shown by the thick solid line. 



meson mass m\a{(3, K) using 



— (/?, K) = \ (6) 



The lines of constant T /Tpc is determined by the ratio of T jmy to Tpc/niy where Tpc/rriY is obtained by T jm^ at 
Kt on the same LCP. We use an interpolation function, Tpc/rriY = A{1 + i3(mps/»Tiv)^)/(l + (^(mps/my)^) with 
A = 0.2253(71), B = -0.933(17) and C = -0.820(39), obtained in Ref. to evaluate Tp^/my for each mps/my. 
The bold dashed line denoted as Kt{Nt = 4) in Fig. [3] represents the pseudo-critical hne T/Tpc = 1. The thin dashed 
lines represent the results for T/Tpc = 0.8, 1.2, 1.4, 1.6, 1.8, 2.0 at Nt = 4. 

We perform finite temperature simulations on a lattice with a temporal extent = 4 and a spatial extent Ns — 16 
along the LCP's for mps/my = 0.65 and 0.80. The standard hybrid Monte Carlo algorithm is employed to generate 
full QCD configurations with two flavors of dynamical quarks. The length of one trajectory is unity and the step size 
of the molecular dynamics is tuned to achieve an acceptance rate greater than 70%. Runs are carried out in the range 
(3 = 1.50-2.40 at thirteen values of T/Tpc ~ 0.82-4.0 for mps/my = 0.65 and twelve values of T/Tpc ^ 0.76-3.0 for 
TTips/'n^v = 0.80. Our simulation parameters and the corresponding temperatures are summarized in Tab.]!] Because 
the determination of the pseudo critical line is more difhcult than the calculation of T/my, the dominant source for 
the error of T/Tpc in TableUis the overall factor Tpc/my. The number of trajectories for each run after thermalization 
is 5000-6000. We measure physical quantities at every 10 trajectories. 



III. HEAVY QUARK FREE ENERGIES 

A free energy of static quarks on the lattice is described by the correlations of the Polyakov loop: il(x) = 
n^=i f^4('''i x) where the t/^(r, x) G SU(3) is the link variable. With an appropriate gauge fixing, one can define 
the free energy in various color channels separately [1^, H^]: the color singlet QQ channel (1), the color octet QQ 
channel (8), the color anti-triplet QQ channel (3*), and the color sextet QQ channel (6), given as follows. 



-F-,{r,T)/T ^ 
-Fa{r,T)/T _ 



i(Trf7t(x)f](y)), (7) 

i (Trf^t (x)Trf7(y)) - ^ (Trf^t (x)f^ (y)) , (8) 

.-P<^ir,T)/T ^ _L(Trl7(x)Trf7(y)) + l(Trl7(x)f](y)), (9) 

-f3.(.,T)/T ^ i(Trl7(x)Trl7(y))-l(TrO(x)17(y)), (10) 
6 6 



where r = |x — y|. We adopt the Coulomb gauge fixing in this study. 
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TABLE I: Simulation parameters for m-ps/' 





K 




Traj. 


1.50 


0.150290 


0.82(3) 


5000 


1.60 


0.150030 


0.86(3) 


5000 


1.70 


0.148086 


0.94(3) 


5000 


1.75 


0.146763 


1.00(4) 


5000 


1.80 


0.145127 


1.07(4) 


5000 


1.85 


0.143502 


1.18(4) 


5000 


1.90 


0.141849 


1.32(5) 


5000 


1.95 


0.140472 


1.48(5) 


5000 


2.00 


0.139411 


1.67(6) 


5000 


2.10 


0.137833 


2.09(7) 


5000 


2.20 


0.136596 


2.59(9) 


5000 


2.30 


0.135492 


3.22(12) 


5000 


2.40 


0.134453 


4.02(15) 


5000 



= 0.65 (left) and mps/mv = 0.80 (right). 
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143749 


0.80(4) 


6000 


1 


70 





142871 


0.84(4) 


6000 


1 


80 





141139 


0.93(5) 


6000 


1 


85 





140070 


0.99(5) 


6000 


1 


90 





138817 


1.08(5) 


6000 


1 


95 





137716 


1.20(6) 


6000 


2 


00 





136931 


1.35(7) 


5000 


2 


10 





135860 


1.69(8) 


5000 


2 


20 





135010 


2.07(10) 


5000 


2 


30 





134194 


2.51(13) 


5000 


2 


40 





133395 


3.01(15) 


5000 




FIG. 4: Simulation results of the normalized free energies scaled by Tpc for color singlet and octet QQ channels (left) and color 
anti-triplet and sextet QQ channels (right) at mps/mv = 0.65. 



Above Tpc, we introduce normalized free energies (Vi, Vs, Ve, Vs*) such that they vanish at large distances. This is 
equivalent to defining the free energies by dividing the right-hand sides of Eqs. ([7))- pUl) by (TrSl)^. 

The normahzed free energies are shown in Fig. 0] for color singlet and octet QQ channels (left) and color anti-triplet 
and sextet QQ channels (right) for m^^/my — 0.65 and T > Tpc- These for mps/my — 0.80 are also shown in 
Fig. [5l Data of the normalized free energies for all temperatures above Tpc are summarized in Appendix |B] From 
these figures, we find that the inter-quark interaction is "attractive" in the color singlet and anti-triplet channels and 
is "repulsive" in the color octet and sextet channels. We also see that, irrespective of the channels, the inter-quark 
interaction becomes rapidly weak at long distances as T increases, as expected from the Debye screening at high 
temperatures. These behaviors are qualitatively similar to the case of quenched QCD in the Lorenz gauge as reported 
in Ref. ^E^. 

To study the screening effects in each color channel more closely, we fit the free energies by the screened Coulomb 
form: 

FM(r,^) = C(M)^2£^e-™°(^)^ (11) 
r 

where acff{T) and moiT) are the effective running coupling and Debye screening mass, respectively. The Casimir 
factor C{M) = {J2a=i ' ^2) a/ for color channel M is explicitly given by 

= ^{8)^1, C(3*) = -| (12) 
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FIG. 5: The same figure with Fig. |4]at mps/my = 0.80. 



TABLE II: The x^/Ndf for each color channel 



T/Tpc 


M = 1 


8 


6 


3* 


1.00 


1.97 


2.26 


2.50 


1.41 


1.07 


2.06 


0.99 


1.04 


1.59 


1.18 


2.08 


1.86 


1.21 


1.59 


1.32 


1.19 


4.02 


1.58 


0.85 


1.48 


1.32 


1.55 


1.06 


0.98 


1.67 


2.96 


1.65 


1.98 


2.64 


2.09 


3.51 


2.16 


2.07 


1.64 


2.59 


1.49 


0.98 


1.25 


3.55 


3.22 


2.05 


1.66 


2.14 


1.83 


4.02 


2.75 


1.51 


1.37 


2.73 



" Fit is unstable since Ve{r,T) at this parameter point is 



temperature at mps/mv = 0.65 (left) and 0.80 (right). 



T/Tpc 


M = 1 


8 


6 


3* 


1.08 


1.25 


1.14 


a 


0.82 


1.20 


1.78 


1.66 


1.33 


0.95 


1.35 


1.32 


0.76 


2.10 


0.81 


1.69 


1.01 


1.58 


1.07 


1.22 


2.07 


4.07 


1.51 


2.28 


2.80 


2.51 


2.00 


2.72 


1.43 


1.93 


3.01 


2.13 


2.60 


1.56 


2.08 



than the statistical errors in the fit range. 



for our cases. Here, it is worth stressing that, with the improved actions we adopt, the rotational symmetry is well 
restored in the heavy quark free energies [Slj . Therefore, we do not need to introduce terms correcting lattice artifacts 
at short distances in Eq. (fTT|) to fit the data shown in Fig. [4] and [5l 

The Debye screening effect is defined through the long distance behavior of Vm {r, T) . A discussion to determine 
an appropriate fit range is given in Appendix VK\ and we choose it to be a/TT/4 < rT < 1.5. We also discuss the 
systematic errors due to a difference of the fit ranges in Appendix \^ and find that the systematic errors are smaller 
than the statistical errors at T ^ 2Tpc- The fit is performed by minimizing x^/^df, where N^p = 20. Results of 
the x^l^DF for each color channel and temperature are summarized in Tab. HIl for mps/my = 0.65 (left) and 0.80 
(right). 

The results of acff{T) and mrtiT) are shown in Fig. [S] for mps/my — 0.65 and Fig. [7] for 0.80. The numerical 
results are also summarized in Tab. Illll for mps/mv = 0.65 and Tab. II VI for 0.80. We find that there is no significant 
channel dependence in acs{T) and m]j(T) at sufficiently high temperatures (T ^ 2Tpc). In other words, the channel 
dependence in the free energy can be well absorbed in the kinematical Casimir factor at high temperatures, as first 
indicated in quenched studies [l3, [ll| . 



A. Debye mass on the lattice and that in perturbative theory 

Let us first compare the Debye mass on the lattice with that calculated in the thermal perturbation theory. The 
2-loop running coupling is given by 

5,f(^)=;3oln(^)%^lnln(i^)', (13) 



a; /3o va 
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FIG. 6: The effective running coupling Qcff(r) (left) and Debye screening mass mD{T) (right) for each color channel as a 
function of temperature from the large distance behavior of the potentials at mps/mv = 0.65. 
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FIG. 7: The same figures with Fig. [6] at mps/my = 0.80 
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TABLE III: Results of acB{T) and mr)(T) at mpg/mv = 0.65 with statistical errors determined by a jackknife method with 
the bin-size of 100 trajectories. 







a 










M = 1 


8 


6 


3* 


1 


8 


6 


3* 


1.00 


2.62(211) 


-0.47(21) 


-0.60(31) 


1.05(67) 


5.00(89) 


2.00(46) 


3.03(56) 


3.63(69) 


1.07 


1.11(43) 


-0.17(6) 


-0.01(0) 


2.98(133) 


4.28(46) 


1.74(44) 


0.09(40) 


5.18(52) 


1.18 


1.49(28) 


0.03(2) 


0.90(50) 


1.38(33) 


4.52(21) 


2.16(79) 


4.75(60) 


4.46(27) 


1.32 


2.03(32) 


0.28(11) 


1.02(27) 


2.24(39) 


4.97(18) 


3.46(43) 


4.52(30) 


5.05(21) 


1.48 


1.12(14) 


0.85(62) 


0.81(26) 


1.14(13) 


4.23(15) 


4.20(83) 


4.04(38) 


4.19(13) 


1.67 


0.69(5) 


0.63(18) 


0.99(20) 


0.72(6) 


3.57(8) 


3.75(32) 


4.19(25) 


3.61(10) 


2.09 


0.83(6) 


0.61(9) 


0.69(6) 


0.66(4) 


3.72(8) 


3.56(17) 


3.65(11) 


3.46(8) 


2.59 


0.71(4) 


0.53(6) 


0.60(5) 


0.72(4) 


3.54(7) 


3.35(14) 


3.44(10) 


3.59(7) 


3.22 


0.55(2) 


0.49(5) 


0.45(2) 


0.48(2) 


3.19(6) 


3.15(11) 


3.05(6) 


3.05(5) 


4.02 


0.48(1) 


0.42(2) 


0.44(1) 


0.51(2) 


3.03(3) 


3.01(8) 


3.02(5) 


3.11(4) 
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TABLE IV: The same Table as Tab. |lll]at mps/mv = 0.80. 







OLcB 


(T) 






mo 


(T) 






M = 1 


8 


6 


3* 


1 


8 


6 


3* 


1.08 


1.23(16) 


-0.03(1) 


— 


1.15(16) 


4.09(16) 


0.63(39) 


— 


3.95(17) 


1.20 


1.61(16) 


2.40(209) 


1.64(48) 


1.59(22) 


4.48(12) 


5.84(102) 


4.87(35) 


4.42(17) 


1.35 


1.18(12) 


0.61(19) 


0.96(17) 


1.13(14) 


4.14(12) 


3.74(36) 


4.10(20) 


4.09(15) 


1.69 


0.78(4) 


0.48(7) 


0.51(5) 


0.81(8) 


3.62(7) 


3.33(16) 


3.33(11) 


3.67(11) 


2.07 


0.64(3) 


0.69(10) 


0.67(6) 


0.63(3) 


3.40(6) 


3.60(17) 


3.51(10) 


3.37(6) 


2.51 


0.56(2) 


0.57(5) 


0.60(4) 


0.55(1) 


3.19(4) 


3.31(11) 


3.36(9) 


3.20(4) 


3.01 


0.51(2) 


0.39(3) 


0.46(2) 


0.53(2) 


3.13(6) 


2.93(10) 


3.05(5) 


3.17(6) 
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FIG. 8: The Debye screening masses moiT) at mps/mv = 0.65 and 0.80 in the color singlet channel together with that 
calculated in the leading-order (left) and next-to-leading-order (right) thermal perturbation theory shown by the dashed lines. 
^ is the renormalization point chosen at = ttT, 2nT, 3ttT. 



where the argument in the logarithms can be written as /i/A = {ii/T){T/Tpc)(Tpc/ A) where we adopt A — ^ ~ 
261 MeV [11] and Tp^ ~ 171 MeV Q. We assume that the renormahzation point fi is in the range fi = ttT to 3nT. 
Therefore, g2i can be viewed as a function of T/Tpc- With g2\, the Debye screening mass is given by 



T 



521 (/^), 



(14) 



in the leading-order (LO) thermal perturbation theory, neglecting the effects of quark masses. 

In Fig. [5] (left) we compare tod(T) in the color singlet channel with m]^ {T)/T for /i = ttT, 2ttT and 3ttT. We find 
that the leading-order screening mass m^jP (T) does not reproduce the simulation results at all. Similar discrepancy 
has been observed also in quenched QCD [s^] and in full QCD with staggered quarks [l2| ■ 

To study higher-order contributions in the thermal perturbation theory, we test the Debye mass in the next-to- 
leading-order calculated by the hard thermal loop resummation method [3J], 



,NLO 
'■D 



(T) 



T 



Nf 



321 (m) 



1 +52i(m) 




^mag 



(15) 



Here mmag(T') — C'mff2i(/^)-^ denotes the magnetic screening mass. The factor Cm cannot be determined within the 
perturbation theory due to the infrared problem. In this study, we adopt C,„ ~ 0.482 obtained by a quenched lattice 



simulation [33| as a typical value ^. In Fig. [8] (right), our simulation results for m/j are compared with 



,NLO 
I'D 



(T) 



^ When we fit Cm from our simulation results m,£){T = 4.02Tpc) with /i = 27rT, we obtain Cm — 0.33. 
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FIG. 9: The ratio R{T) in the text which is supposed to be close to unity if moiT) = yjl + Nf/6 geff (T) T holds. The results 
are for the color singlet channel. 



shown by the dashed lines. The values of m^'"*-* are larger than m]^ by approximately 50 % and agree better with 
the simulation results. 



B. Phenomenological relation between acff and mo 

So far, we have fitted the free energies on the lattice treating acfr and mjj as independent parameters. In this 
subsection, we test an ansatz inspired from the leading-order perturbation theory: 



,16) 



with gcs(T) = 47racff (T). Therefore, if this relation holds, we expect that the ratio 

RiT) . ^^ + ^fl^)-'''-^o^T)/T ^^^^ 
A/47raoff(r) 

should be close to unity [R{T) ^ 1). 

In Fig. [HI our simulation results of R{T) for the singlet channel are shown as a function of T/Tpc- We find that R{T) 
is consistent with unity for T ^ l.bTpc within 10% accuracy. This is a non-trivial observation particularly near Tpc 
and suggests that the major part of the higher-order effects and non-perturbative effects of moiT) can be expressed 
by the effective running coupling gcs{T). We note that a similar effective coupling at T = defined through the 
lattice potential was introduced to improve the lattice perturbation theory • 



C. Comparison with the staggered quark action 

Finally, we compare the results of aoff(T) and m£)(r) obtained by the Wilson quark action (present work) with 
those by an improved staggered quark action on a 16^ X 4 lattice at mps/my =i 0.70 [Hlil- The comparison is shown 
in Fig. [To] for acff(r) (left panel) and mo{T) (right panel). Although a^siT) does not show significant difference 
between the two actions, mz){T) in the Wilson quark action is systematically higher than that of the staggered quark 
action by about 20% even at T = 4Tpc- This discrepancy can be seen directly from the normalized free energies in Fig. 
[TT] where —r x Vi(r) is shown as a function of rT at T ~ l.lTpc (left) and ~ 4.0Tpc (right). The circles (triangles) are 
the results of the Wilson quark action (staggered quark action). The data for the staggered quark action is taken from 
Ref. The straight lines are the fits with the screened Coulomb form, Eq. ([TT|) . in the range vTT/4 < rT < 1.5 
(circles) and 0.8 — 1.0 ^ rT (triangles). The intercepts of the lines with the vertical axis and the slopes of the lines 
correspond to aeff(T) and mjjiT), respectively. There is an obvious difference in the slope in different quark actions, 
which may be regarded as systematic errors due to the lattice discretization. This discrepancy should be further 
investigated at smaller lattice spacing, i.e. larger lattice size in the temporal direction. 
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FIG. 10; Comparison of the aeff (r)(left) and m_D(r) (right) between the results of the Wilson quark action and staggered quark 
action. 




FIG. 11: Comparison of the normalized free energies, — r x Vi(r), between the Wilson quark action (present work) and the 
staggered quark action in Ref. [U as a function of rT at T ~ 1.7Tpc (left) and ~ 4.0Tpc (right) with the log scale. Straight 
lines show the fits based on the screened Coulomb form in the range \/TT/4 < rT < 1.5 (circles) and 0.8 — 1.0 ^ rT (triangles). 



D. Force between heavy quarks 



To make a direct comparison of the free energies below and above Tpc, we study the "force" between heavy quarks 
defined by dFM{r,T)/dr without introducing the subtraction of (Trl7)^. In Fig. [121 such forces are shown for color 
singlet (left) and octet (right) QQ channels, and in Fig. [131 for color sextet (left) and anti-triplet (right) QQ channels 
at mps/m-v = 0.65. Results at mps/my — 0-80 are shown in Fig. [Hland[T5l 

In the color singlet and anti-triplet channels, there is always an "attraction" both for T < Tpc and for T > Tpc, and 
the attraction is stronger as we decrease the temperature. The signal for the color sextet channel below Tp^ is rather 
weak and we cannot make a definite statement from our data. On the other hand, the force in the color octet channel 
is clearly "attractive" ("repulsive") at low (high) T. Namely, the simple Casimir scaling law, which worked well at 
high T, does not hold below Tpc- In Ref. [3y|, it has be suggested that the Polyakov loop correlation in the octet 
channel is saturated by the color-singlet intermediate states at very low temperature. If this is true, the attraction in 
the color octet channel found here is explained by the attraction in the singlet channel. 



IV. SPATIAL WILSON LOOP 

In previous quenched studies [13, HI, [I^l , the Wilson loop in spatial direction is found to show non- vanishing 
spatial string tension as even at T > Tpc, which is called the spatial confinement. We study this phenomenon with the 
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FIG. 12: Results of dFM{r)/dr scaled by Tpc for color singlet (left) and octet (right) QQ channels at mps/rny = 0.65. 




FIG. 13: Results of dFAr{r)/dr scaled by Tpc for color sextet (left) and anti-triplet (right) QQ channels at mps/mv ~ 0.65. 



existence of dynamical quarks in the system. Although quarks are expected to decouple from the spatial observables 
for T 3> Tpc due to dimensional reduction and thus do not affect cts in the high temperature hmit, it is not obvious 
whether the same is true near Tpc. 

We evaluate as assuming a simplest ansatz for the spatial Wilson loop W{I, J) with the size I x J: 

-\nW{I,J)^aJJ + crp{2I + 2J) + C^, (18) 

where as, ap and Cw are fit parameters. The results of -y/ as (T) /Tpc are shown in Fig. [16] (left) as a function of T/Tpc. 
We find that \/ajT)/Tpc approaches to a constant as T decreases below Tpc, while it increases linearly as T above 
Tpc- Similar behavior has been observed in the quenched case [13] ■ 

Let us first make a phenomenological parameterization of the spatial string tension: 

y^:(T) = cgl,{fi)T, (19) 

motivated by the dimensional reduction at high temperature with g|i(/i) being the two- loop running coupling of two- 
fiavor QCD defined in Eq. p^ . Our results for ^Jas{T)/T are shown in Fig. [TBI (right). The solid and dashed lines 
in the figure are the results of a fit of our data in the range T > l.STpc to the formula Eq. p^ . Here, we fix /i = 2TrT 
and take c and A/Tpc as fitting parameters to obtain 



c = 0.694(12), A/Tpc = 0.866(48) for mps/my = 0.65, 
c = 0.639(34), A/Tpc = 0.92(14) for mps/my = 0.80. 



(20) 
(21) 
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FIG. 14: The same as Fig.[T2]at mps/my = 0.80. 
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FIG. 15: The same as Fig. [13] at mps/my = 0.80. 



Similar numbers have been observed in a quenched simulation [20], (c, A/Tpc) — (0.566(13), 0.653(57)), and in a 2 + 1 
flavor simulation with staggered quark [37], (c, A/Tpc) = (0.587(41), 0.72(17)). (Here we rescaled the previous results 
of A by 2tt to compare at the common value oi fj, = 2'itT .) 

Let us now make an alternative comparison of our data with a recent prediction by the parameter-free three- 
dimensional (3d) eS'ective theory [2l|, IsB] , which gives 

= 0.553(1) 5^, (22) 

where the coefficient 0.553(1) expresses a non-perturbative contribution determined by Zd quenched lattice simulations, 
and gyi is a dimensionful id gauge coupling defined through the four-dimensional running coupling, g{^)^ as 

9m = .921 (m) T[l + C^ gl ifi) + C2 gt, (m) + ' ' ' ] • (23) 

Known coefficients Ci and C2 represent the next-to-leading (NLO) and next-to-next-leading (NNLO) contributions 
in 3d effective theory [2l|. We choose fi/T — 2tt, which is in the range 6.0-7.0 assumed in the study of Ref. [IJ. We 

take the same value for A/Tpc — ^j^'^ f^pc as that in Section ITlI Al 

Figure [TTl shows the results of the spatial string tension, together with the prediction from the 3c? effective theory 
in the leading order, NLO and NNLO. We observe that convergence of the series is rather well and the NNLO result 
is quite consistent with our lattice data. Similar agreement in quenched QCD has been reported in Ref. ,21i] . 
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FIG. 16: The square root of the spatial string tension over Tpc (left) and T (right) as a function of T/Tpc at mps/mv ~ 0.65 
and 0.80. The solid and dashed lines in the right panel express the fit results based on Eq. (|19|l at mps/mv = 0.65 and 0.80, 
respectively. 




FIG. 17: Comparison of our results and predictions of the three-dimensional effective theory with contributions up to leading 
order (solid), NLO (dashed) and NNLO (dash-dotted line), respectively. 



V. CONCLUSIONS 



In the past few years, many of the lattice QCD simulations at finite temperature and density have been performed 
using staggered quark actions. However, to control the lattice artifacts, comparison with other quark actions are 
indispensable. This motivates us to carry out a systematic study with an improved Wilson quark action. As a first 
step, we performed simulations of Nf — 2 QCD on an x Nt ^ 16^ x 4 lattice. We have identified the lines of 
constant physics and studied the temperature-dependence of various quantities at mps/my = 0.65 and 0.80 in the 
range T/Tpc - 0.76-4.0. 

We found that, at T > Tpc, the free energies of QQ and QQ normalized to be zero at large separation show 
attraction (repulsion) in the color singlet and anti-triplet channels (color octet and sextet channels). We fitted the 
free energy data in each channel by the screened Coulomb form, Eq. pT|) . which consists of the Casimir factor, the 
effective coupling (T) and the Debye screening mass m d (T) . We found that free energies in different channels can 
be fitted for T ^ 2Tpc in terms of universal acs{T) and m£){T) while all the channel dependence can be absorbed 
by the Casimir factor. The magnitude and the T-dependence of m£)(T) is consistent with the next-to-leading order 
calculation in thermal perturbation theory. Moreover it is also well approximated by the leading order form with an 
"effective" running coupling defined from acs{T). 

By comparing our results with the improved Wilson quark action and those with the improved staggered quark 
action, we found that acii{T) does not show appreciable difference while ttidIT) in the Wilson quark action is larger 
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than that of the staggered quark action by 20%. To draw a definite conclusion, however, simulations with smaller 
lattice spacings, i.e., larger lattice sizes in the temporal direction (such as Nt — 6 or larger) at smaller quark masses 
are required. 

We also discussed the force between heavy quarks to make a direct comparison of the free energies below and above 
Tpc- In the color octet channel, we find that the force becomes attractive below Tpc against the simple Casimir scaling 
law. This may be related to the suggestion that the Polyakov loop correlation in the octet channel is not independent 
from that in the singlet channel in the low temperature limit [301. Moreover, spatial string tension <Js{T) in the 
high temperature phase was investigated. The result of the spatial string tension in the quark-gluon plasma shows 
a behavior consistent with a/ct7(T) = c(72i(2vrT)r and agrees well with a parameter-free prediction of the NNLO 
three-dimensional effective theory. 

Using our configurations, we are currently studying the critical temperature (T^j, in the chiral limit), the chiral and 
isospin susceptibility across the phase transition, the effect of the chemical potential on the heavy-quark free energies, 
and the equation of state at finite temperature and density. Some of the results are reported in Ref. (39l . |40| . 
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APPENDIX A: FIT RANGE AND SYSTEMATIC ERRORS 



In order to determine the appropriate fit range for free energies by the screened Coulomb form, Eq. (jlip , we estimate 
the effective Debye mass from the ratio of normalized free energies: 



Ar) 



■log 



Ar 

r 



(Al) 



The effective Debye mass is shown in Fig. [T51 for color singlet channel at T = 1.18Tpc (left) and 3.22Tpc (right) for 
mps/my = 0.65. Since the plateaus of rrf^{T; r) are found at rT > 0.8, we choose the fit range to be \/TT/4 < rT < 
1.5. The solid lines in Fig. [18] are the results of the fit with the statistical errors indicated by the dashed lines. 

In order to investigate systematic errors due to a dependence on fit ranges, we consider two other fit ranges: 
V9/4 < rT < 1.5 and Vl3/A < rT < 1.5. The systematic errors due to the difference of the lower end of the fit range 
are shown in Fig.[Tn]for acff(T) (left) and m£)(T) (right) in the color singlet channel at mps/my = 0.65. The similar 
plots at mps/mv = 0.80 are shown in Fig. [201 We find that, for T ^ 2Tpc, the systematic errors amount to at most 
the twice of the statistical errors, i.e. 2a. On the other hand, for T > 2Tpc, the systematic errors are less than 10% 
of the statistical errors. 
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FIG. 18: The effective Debye mass defined by a ratio of normalized free energies, Eq. (|A1|| . for color singlet channel at 
T = l.lSTpc (left) and 3.22rpc (right) for mps/mv = 0.65. The fit results (statistical errors) are also given as solid lines 
(dashed lines). 
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FIG. 19: The effective running coupling (left) and Debye screening mass (right) for color singlet channel with the systematic 
errors due to the difference of the fit range defined in text at mps/mv = 0.65. 
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APPENDIX B: DATA LISTS OF NORMALIZED FREE ENERGIES 

In this Appendix, we give data of the normalized free energies for aU color channels and temperatures above Tpc 
as a function of distance (rT). The normalized free energies are measured at every ten trajectories for each quark 
mass and temperature summarized in Tab. IH The values in parenthesis of the normalized free energies express the 
statistical errors determined by a jackknife method with the bin-size of 100 trajectories. 



TABLE V: Data lists of normalized free energies for all color channels at T = l.OOTpc (left) and 1.07Tpc (right) for mps/my = 
0.65 as a function of rT 





r = l.OOTpc, mps/mv = 0.65 


T = 1.07Tpe, mps/mv = 0.65 


rT 


Vi 




Ve 


Vs- 


Vi 


Vs 


Ve 


Vs- 




0.250 


-2.3976(132) 


-0.1786(59) 


0.0466(78) 


-1.0760(69) 


-1.6513(71) 


0.0337(27) 


0.1796(31) 


-0.7895 


'34) 


0.354 


-1.2753(110) 


-0.1199(51) 


-0.0228(56) 


-0.5567(58) 


-0.8627(51) 


0.0041(18) 


0.0783(24) 


-0.4158 


;28) 


0.433 


-0.7376(104) 


-0.0931(60) 


-0.0390(75) 


-0.3389(72) 


-0.5205(46) 


-0.0013(22) 


0.0430(28) 


-0.2525 


;27) 


0.500 


-0.5979(95) 


-0.0956(74) 


-0.0471(90) 


-0.2766(90) 


-0.4137(52) 


-0.0085(25) 


0.0251(26) 


-0.2042 


;28) 


0.559 


-0.3640(72) 


-0.0618(39) 


-0.0352(42) 


-0.1726(54) 


-0.2728(33) 


-0.0077(11) 


0.0152(13) 


-0.1388 


19) 


0.612 


-0.2347(64) 


-0.0503(38) 


-0.0384(36) 


-0.1287(54) 


-0.1884(32) 


-0.0087(13) 


0.0076(15) 


-0.0997 


21) 


0.707 


-0.1170(118) 


-0.0449(54) 


-0.0405(56) 


-0.0824(60) 


-0.1117(39) 


-0.0108(14) 


0.0009(20) 


-0.0638 


'21) 


0.750 


-0.1010(77) 


-0.0318(42) 


-0.0276(42) 


-0.0577(51) 


-0.0860(29) 


-0.0091(11) 


-0.0015(12) 


-0.0507 


16) 


0.791 


-0.0727(82) 


-0.0241(34) 


-0.0220(40) 


-0.0425(50) 


-0.0727(36) 


-0.0073(15) 


-0.0016(17) 


-0.0418 


'18) 


0.829 


-0.0566(109) 


-0.0223(35) 


-0.0098(45) 


-0.0405(53) 


-0.0553(40) 


-0.0072(12) 


-0.0020(13) 


-0.0324 


;i7) 


0.866 


-0.0719(161) 


-0.0232(55) 


-0.0169(57) 


-0.0485(101) 


-0.0332(56) 


-0.0037(22) 


0.0006(25) 


-0.0279 


;3o) 


0.901 


-0.0364(97) 


-0.0140(46) 


-0.0170(41) 


-0.0300(60) 


-0.0361(29) 


-0.0062(10) 


-0.0028(12) 


-0.0214 


'20) 


0.935 


-0.0401(72) 


-0.0136(31) 


-0.0171(33) 


-0.0269(38) 


-0.0295(26) 


-0.0059(8) 


-0.0033(11) 


-0.0191 


'14) 


1.000 


-0.0045(212) 


-0.0151(74) 


-0.0007(92) 


-0.0046(114) 


-0.0234(54) 


-0.0042(22) 


-0.0008(23) 


-0.0152 


;35) 


1.031 


-0.0157(64) 


-0.0138(26) 


-0.0115(32) 


-0.0184(39) 


-0.0147(24) 


-0.0040(11) 


-0.0016(11) 


-0.0108 


;i5) 


1.061 


-0.0203(71) 


-0.0103(27) 


-0.0093(32) 


-0.0172(52) 


-0.0177(28) 


-0.0038(9) 


-0.0038(12) 


-0.0112 


;i8) 


1.090 


-0.0305(92) 


-0.0104(36) 


-0.0053(40) 


-0.0063(67) 


-0.0082(38) 


-0.0033(13) 


-0.0027(13) 


-0.0066 


'18) 


1.118 


0.0008(90) 


-0.0158(35) 


-0.0140(50) 


-0.0206(58) 


-0.0112(32) 


-0.0048(11) 


-0.0025(13) 


-0.0051 


'17) 


1.146 


-0.0103(74) 


-0.0073(25) 


-0.0049(33) 


-0.0097(38) 


-0.0048(19) 


-0.0032(10) 


-0.0025(9) 


-0.0057 


'14) 


1.173 


-0.0149(125) 


-0.0109(45) 


-0.0027(36) 


-0.0027(64) 


-0.0064(31) 


-0.0034(13) 


-0.0022(13) 


-0.0072 


;i8) 


1.225 


-0.0024(101) 


-0.0056(31) 


-0.0100(39) 


-0.0127(49) 


-0.0075(31) 


-0.0029(11) 


-0.0027(15) 


-0.0066 


;i8) 


1.250 


-0.0136(82) 


-0.0068(33) 


-0.0036(33) 


-0.0127(57) 


-0.0065(29) 


-0.0028(11) 


-0.0016(12) 


-0.0045 


;i7) 


1.275 


-0.0084(61) 


-0.0047(18) 


-0.0043(20) 


-0.0080(33) 


-0.0082(17) 


-0.0018(8) 


-0.0011(9) 


-0.0044 


;i2) 


1.299 


0.0162(83) 


-0.0041(35) 


-0.0043(35) 


-0.0050(52) 


-0.0045(24) 


-0.0028(9) 


-0.0022(11) 


-0.0059 


'15) 


1.346 


-0.0008(63) 


-0.0038(22) 


-0.0010(23) 


-0.0038(31) 


-0.0043(19) 


-0.0019(6) 


-0.0011(7) 


-0.0056 


'12) 


1.369 


-0.0070(79) 


-0.0074(26) 


-0.0075(31) 


-0.0062(39) 


-0.0038(19) 


-0.0006(8) 


-0.0013(9) 


-0.0034 


'13) 


1.414 


-0.0106(159) 


-0.0150(44) 


0.0006(47) 


-0.0137(88) 


-0.0021(53) 


-0.0015(13) 


-0.0001(16) 


-0.0037 


'26) 


1.436 


0.0013(73) 


-0.0029(27) 


-0.0006(34) 


-0.0078(41) 


-0.0042(22) 


-0.0008(8) 


-0.0005(10) 


-0.0019 


14) 


1.458 


-0.0112(70) 


0.0013(28) 


-0.0010(31) 


-0.0095(47) 


-0.0015(24) 


-0.0003(10) 


-0.0004(11) 


-0.0014 


.15) 
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-0.0034(70) 


-0.0036(29) 


-0.0014(29) 


0.0015(45) 


-0.0012(22) 


-0.0010(8) 


-0.0010(9) 


-0.0015 


'12) 


1.500 


0.0178(119) 


0.0020(31) 


0.0011(35) 


-0.0078(62) 


-0.0048(24) 


-0.0013(9) 


-0.0006(12) 


-0.0031 


'16) 
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TABLE VI: Data lists of normalized free energies for all color channels at T = l.lSTpc (left) and 1.32Tpc (right) for mps/mv = 
0.65 as a function of rT 
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1.414 
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0.0009(16) 
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0.0001(7) 
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0.0007(12) - 
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-0.0003(4) 
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1.458 
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-0.0015(10) 


0.0008(4) 


0.0006(4) 


-0.0010(5) 
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-0.0037(12) 


-0.0003(5) 


-0.0007(5) 


-0.0008(6) 


-0.0011(11) 


0.0006(4) 


0.0006(4) 


-0.0006(5) 


1.500 


-0.0040(19) 


-0.0009(6) 


-0.0012(7) 


0.0001(9) 


0.0005(12) 


0.0004(3) 


0.0009(4) 


0.0001(8) 
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TABLE VII: Data lists of normalized free energies for all color channels at T = 1.48Tpc (left) and 1.67Tpc (right) for mps/mv = 
0.65 as a function of rT 





T 
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T 
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u. i iu / yyj ^ 


U.zoDDl^iz J 


n r^99(' 1 9^ 
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TABLE VIII: Data lists of normalized free energies for all color channels at T = 2.09Tpc (left) and 2.59Tpc (right) for mps/ 
0.65 as a function of rT 
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TABLE IX: Data lists of normalized free energies for all color channels at T = 3.22Tpc (left) and 4.02Tpc (right) for mps/ 
0.65 as a function of rT 
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TABLE X: Data lists of normalized free energies for all color channels at T = l.OSTpc (left) and 1.20Tpc (right) for mps/mv = 
0.80 as a function of rT 
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TABLE XI: Data lists of normalized free energies for all color channels at T = 1.35Tpc (left) and 1.69Tpc (right) for mps/mv = 
0.80 as a function of rT 
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TABLE XII: Data lists of normalized free energies for all color channels at T = 2.07Tpc (left) and 2.51Tpc (right) for mps/ 
0.80 as a function of rT 
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TABLE XIII: Data lists of normalized free energies for all color channels at T = S.OlTpc for mps/mv = 0.80 as a function of 
rT 
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